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The  Faraday  rotation  and  amplitude  data  for  the  period  September  1978  through  August  1979, 
recorded  from  the  geostationary  satellite  SMS  I  |90°W)  at  Natal  (5.9°S,  3,S,2‘’W),  have  been  utilized  to 
study  the  behavior  of  total  electron  content  (TEC).  Apart  from  diurnal  and  seasonal  behavior  of  TEC, 
two  characteristic  features  of  the  nighttime  ionosphere  arc  discussed:  1 1 )  a  posisunset  enhancement  in 
TFiC  which  lasts  for  several  hours  following  the  rapid  sunset  decrease  and  (2)  the  sharp,  isolated 
changes  in  TEC.  In  particular,  relating  to  phenomenon  2.  the  depletions  in  TEC  are  usually 
accompanied  by  a  simultaneous  increase  in  fading  rate,  scintillation  index,  and  amplitude  level,  while 
the  enhancements  are  accompanied  by  a  simultaneous  decrease  in  fading  rate,  scintillation  index,  and 
amplitude  level.  The  statistics  of  their  cKcurrence  and  nature  are  described.  The  average  behavior  of 
the  signal  amplitude  after  propagating  through  cither  a  depleted  or  an  enhanced  ionosphere  is  modeled 
theoretically  and  shown  to  be  in  agreement  with  the  observed  experimental  behavior. 


Behavior  of  Total  Electron  Content 
(TEC)  AT  Natal 


Diurnal  and  Sea,s'onal  Variation  in  TEC 


Natal  is  an  equatorial  station  (latitude  5.9°S,  dip  -9.6°).  A 
rotation  of  tt  rad  on  the  136.379-MHz  signal  received  from 
the  geostationary  satellite  SMS  I  (90°W)  is  equivalent  to  1 .89 
•X  lO”  el/m’  in  total  electron  content  for  a  mean  field  height 
of  350  km.  To  study  the  general  behavior  of  TEC  in  different 
seasons,  10-12  days  of  uninterrupted  data  were  selected 
from  the  periods  September-October  1978,  December  1978 
to  January  1979,  March-April  1979,  and  June-July  1979.  The 
hourly  median  values  of  TEC  are  plotted  in  Figure  I.  it  may 
be  pointed  out  here  that  no  claims  are  made  regarding  the 
absolute  value  of  the  TEC,  especially  for  individual  cases 
near  the  diurnal  minimum  values,  for  want  f  independent 
calibration.  It  is,  however,  felt  that  the  error  in  the  median 
values  is  only  a  small  fraction  of  tr  and  that  no  nw  ambiguity 
is  involved,  though  it  may  sometimes  be  present  in  indivdual 
day-to-day  values. 

Diumally,  the  minimum  occur  around  05(X)  hours  LMT  in 
all  the  seasons,  while  the  maximum  occurs  around  1500 
LMT  in  the  June  solstice  (winter),  1200  LMT  in  the  Decem¬ 
ber  solstice  (summer),  1330  LMT  in  the  spring  equinox,  and 
1430  LMT  in  the  fall  equinox.  The  electron  content  increases 
rapidly  following  sunrise.  The  rate  of  increase  is  very  fast 
during  spring,  fall,  and  the  December  solstice  but  very  slow 
during  the  June  solstice.  After  attaining  the  diurnal  maxi¬ 
mum  value  the  electron  content  decreases  very  fast  during 
spring,  fall,  and  the  June  solstice  but  rather  slowly  during  the 
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December  .solstice.  The  average  diurnal  minimum  and  maxi 
mum  values  of  TEC  for  spring,  fall,  winter,  and  summer  are 
1.33  X  10”  and  9.07  x  I0”.  1.11  x  lO”  and  11.47  x  |0”, 
0.67  X  10”  and  5.70  x  10”,  and  1.00  x  10”  and  7.94  x  lO” 
m  ■  respectively,  giving  a  diurnal  ratio  of  TEC  of  6.8.  10.5, 
8.8,  and  7.9.  It  may  be  noted  that,  in  general,  the  decay  rate 
of  the  TEC  after  sunset  is  greater  than  the  buildup  rate  at 
sunrise. 

Po.yl.yunsel  Enhancement  in  TEC 

Following  the  dramatic  drop  in  Faraday  rotation  angle  for 
an  hour  or  so  after  sunset,  it  may  then  either  settle  down  to  a 
gradual  decay,  finally  reaching  a  minimum  just  before  the 
following  sunrise,  or  sometimes,  and  at  different  seasons,  a 
substantial  increase  in  TEC  occurs  for  a  duration  of  3-4 
hours  centered  around  23(X)  LMT.  This  phenomenon  is 
much  more  pronounced  during  equinoxes  and  hence  is  very 
visible  even  in  median  curves,  while  only  a  trace  increase  is 
noticeable  in  solstices.  On  the  evening  of  November  25, 
1978.  as  shown  in  Figure  2,  the  increase  at  2200  LMT  is 
approximately  L657r  rad  or  3.12  x  10”  el/m‘,  which  is 
nearly  one  half  of  the  diurnal  variation,  while  on  the  night  of 
November  15,  1978,  the  increase  at  2230  LMT  is  3.8w 
radians  or  7.18  x  lO”  el/m’,  which  is  about  90%  of  the 
diurnal  maximum  and  about  4%  more  than  the  diurnal 
change  >n  TEC  (Figure  3).  This  nocturnal  increase,  in 
varying  degrees,  has  been  observed  both  on  days  of  no 
scintillation  as  well  as  on  days  of  strong  scintillation.  There¬ 
fore  it  is  not  likely  to  be  connected  with  the  irregularity- 
producing  mechanism. 

The  large  postsunset  enhancements  in  equatorial  zone 
have  been  reported  by  various  workers  (e.g..  Hunter,  1%9; 
Koster,  1972;  Koster  and  Beer,  1972;  Yehoah-Amankwah 
and  Koster,  1972;  Donatelli  et  al,,  I981J.  It  was  suggested 
that  the  F  region  dynamo  theory  [Rishbeth,  1971,  1977)  may 
offer  the  most  plausible  explanation  for  this  effect.  Rishbeth 
1 1971,  1977]  suggested  that  the  neutral  winds  in  the  F  region 
produce  polarization  fields  that  short-circuit  through  the  E 
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F-'ig.  I.  Diurnal  and  seasonal  variation  of  TEC  over  Natal,  Brazil, 
for  the  period  September  1978  through  August  1979. 


region  in  daytime  but  build  up  at  sunset  through  the  sudden 
drop  in  E  region  conductivity.  This  produces  large  upward 
drifts,  resulting  in  a  rapid  lifting  of  the  F  layer  and  massive 
transport  of  plasma  along  field  lines.  The  rapid  lifting  causes 
an  artihcial  apparent  decrease  in  the  total  Faraday  rotation 
because  the  layer  is  carried  to  higher  altitudes  at  which  the 
magnetic  field  is  weaker,  and  in  addition,  the  upward  trans¬ 
port  causes  a  real  physical  decrease  in  the  vertical  columnar 
TEC  at  the  equator  over  and  above  that  expected  from  the 
artificial  decrease  alone.  The  combined  effects  of  these 
mechanisms  result  in  a  rapid  sunset  decline  in  Faraday 
rotation  as  noted  by  Koster  [1 972)  with  a  dip  latitude  of  4°. 
This  effect  is  quite  pronounced  in  the  present  observations 
from  Natal,  less  pronounced  in  Wun/er’s  11%91  observations 
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Fig.  2.  Typical  Faraday  rotation  variation  showing  postsunset 
enhancement  as  observed  on  November  25,  1978.  This  day  is 
marked  by  the  absence  of  scintillation  even  though  scintillations  are 
usually  very  active  in  November. 
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Fig.  -V  Typical  Faraday  rotation  variation  showing  sharp,  iso¬ 
lated  depletions  superposed  over  the  enhancements  in  the  period 
2100-2330  as  well  as  the  normal  decaying  nighttime  ionosphere  after 
local  midnight. 


at  13°  dip  latitude,  and  even  less  in  the  Ascension  Island 
observations  \Donaielli  el  al.,  1981 )  at  1.5°  dip  latitude.  These 
observations  are  consistent  with  the  Defense  Meteorological 
Satellite  Program  measurements  in  the  postsunset  ( 1930  LT) 
sector  at  840  km  which  showed  smooth  density  depletions  in 
the  immediate  vicinity  of  the  magnetic  equator  with  en¬ 
hancements  in  the  ±l0°-20°  magnetic  latitude  range.  The 
subsequent  increase  around  local  midnight  in  TEC.  howev¬ 
er,  is  more  pronounced  in  the  present  and  in  Hunter's 
observations. 

We  note  that  the  above  arguments  [cf.  Koster.  1972) 
clearly  explain  the  dramatic  decrease  in  TEC  after  sunset  at 
the  magnetic  equator  and  the  relative  increase  in  electron 
density  in  the  10°-20°  magnetic  latitude  range,  but  the  very 
large  enhancements  in  TEC  around  2200-2300  LT  in  the 
equatorial  ionosphere  as  reported  in  this  paper  are  not 
explicitly  explained.  Lanzerotti  et  al.  [1976]  reported,  dis¬ 
cussed,  and  related  the  occurrence  of  nighttime  ledges  in 
TEC  measurements  at  Arecibo,  via  a  model  calculation,  to 
the  rise  and  fall  of  the  height  of  the  Ft  peak  around  midnight. 
They  concluded  that  the  characteristics  of  the  nighttime 
ledge  observed  at  Arecibo  could  be  directly  related  to  the 
meridional  component  of  nighttime  thermospheric  winds. 
Anderson  and  Roble  [1974]  studied  the  behavior  of  the 
nighttime  ionosphere  near  the  geomagnetic  equator  and 
showed  that  the  E  x  B  vertical  drifts  can  cause  perturbations 
in  the  neutral  winds  and  temperatures.  Recently,  Bitiencourl 
and  Sahai  [1979],  using  a  tropical  ionospheric  F  recion 
model  described  by  Bitiencourt  and  Tinslev  [1976],  investi¬ 
gated  the  effect  of  E  x  B  drift  velocity  on  the  local  time 
dependence  of  N^Fi,  h^Fi,  and  6300-A  emission  at  the 
magnetic  equator.  They  concluded  that  the  times  of  N„F2 
enhancements  during  the  night  are  closely  related  to  the  time 
history  of  the  vertical  E  x  B  drift,  its  magnitude,  and  time  of 
reversal.  When  the  postsunset  enhancement  in  upward  drift 
is  well  developed  and  downward  drift  commences  after  2000 
LT,  then  N^Fi  enhancement  occurs  around  0200  LT.  When 
upward  drift  is  less  pronounced  and  downward  drift  occurs 
at  1900  LT,  the  enhancement  occurs  around  2230  LT  [cf. 
Bitiencourt  and  Sahai,  1979;  Anderson  and  Rusch,  1980]. 
This  clearly  explains  the  variability  of  TEC  enhancements 
during  the  nighttime  at  equatorial  stations. 
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Fig.  4.  Original  record  for  November  16,  1978,  from  Natal  (5.9°S,  35,2°W)  using  1.36-MHz  transmissions  from 
satellite  SMS  I  I90'’W)  showing  the  occurrence  of  an  isolated  bubble  and  its  accompanying  increases  in  the  scintillation 
index,  the  average  amplitude,  and  the  fading  rate. 


Isolated  Depletions  and  Enhancements  in  TEC 

The  recent  measurements  at  the  magnetic  equator  have 
helped  to  reach  some  understanding  about  the  nature  of  one 
kind  of  ionospheric  perturbation  occurring  at  night.  Such 
perturbations  are  described  variously  as  ionization  deple¬ 
tions,  bubbles,  plumes,  and  some  type  of  irregularity  patch¬ 
es,  or  equatorial  spread  F  irregularities.  On  the  experimental 
side  the  radar  backscatter  maps  [Woodman  and  LaHoz, 
1976]  show  dramatically  the  plumelike  structures  when  they 
are  viewed  vertically  upward  as  a  function  of  time  at  a  fixed 
location.  In  in  situ  measurements,  large  depletions  of  ioniza¬ 
tion  are  observed  as  the  satellite  passes  nearly  horizontally 
[McClure  et  al.,  1977]  or  as  the  rocket  traverses  nearly 
vertically  [Kelley  et  al.,  1976]  through  these  bubbles.  On  the 
basis  of  propagation  measurements  the  occurrence  of  these 


bubbles  has  been  shown  to  be  correlated  with  the  onset  of 
scintillation  on  satellite  radio  signals  [Aarons,  1977;  Basil 
and  Basu,  1976],  These  have  also  been  found  to  be  responsi¬ 
ble  for  backscattering  in  transequatorial  propagation 
[Rottger,  1976].  On  the  theoretical  side,  early  investigations 
have  shown  that  the  bottomside  ionosphere  is  subject  to 
Rayleigh-Taylor  instabilities  [Liu  and  Yeh,  1966a,  h[.  Recent 
numerical  simulations  have  shown  that  as  the  depletion 
increases,  the  nonlinear  polarization  E  x  B  forces  in  the 
equatorial  geometry  can  drive  the  depleted  region  upward 
through  the  ionization  peak  to  the  topside  ionosphere  |0.y.va- 
kow  et  al.,  1979], 

Yeh  et  al.  [1979a,  h]  reported  observations  of  ionospheric 
bubbles  utilizing  the  136-MHz  Faraday  rotation  records  from 
Natal,  and  Wernik  [1979]  presented  a  simple  model  of  an 
equatorial  bubble.  Looking  again  at  Figure  3,  it  may  be  noted 
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Fig.  6.  Original  record  for  Seplember  19,  1978.  from  Natal  using  SMS  I  transmissions,  showing  an  isolated 

enhancement. 


that  there  are  several  sharp  structures  superposed  over  the 
background  TEC  variation.  The  first  sharp  depletion  oc¬ 
curred  at  1920  LMT,  shortly  after  the  onset  of  scintillations. 
There  are  also  other  depletions,  but  those  occurring  before 
local  midnight  are  not  as  isolated  and  distinct  as  those 
occurring  near  or  after  midnight.  For  example,  the  bubble 
occurring  at  0040  LMT  on  November  16,  1978,  appeared  in 
the  original  record  as  shown  in  Figure  4.  Note  the  coincident 
increase  in  scintillation  index  as  well  as  scintillation  rate  and 
the  average  amplitude  level  with  the  appearance  of  the 
bubble.  This  suggests  that  the  isolated  lapid  scintillation  rate 
can  be  used  as  a  possible  indicator  of  a  bubble.  It  may  be 
mentioned  here  that  these  isolated  bubbles  should  not  be 
taken  as  traveling  ionospheric  disturbances  (TID),  which 
generally  have  quasi-periodic  and  wavelike  structure.  It  is 
very  clear  from  Figure  3  that  the  five  sharp  depletions 
occurring  successively  a  few  hours  around  local  midnight  do 
not  possess  the  usual  TID  appearance. 

Another  example  is  shown  in  Figure  5.  The  amplitude 
channel  at  the  bottom  shows  rapid  fluctuations  at  three 
times:  near  0015,  0025,  and  0050  LMT.  Coincidentally,  the 
depletions  in  Faraday  rotation  angle  can  be  clearly  identified 
for  the  latter  two  times  but  not  so  clearly  for  the  former  time. 
Every  time  the  average  amplitude  level  also  shows  an 
increase.  In  addition  to  sudden  increases  in  scintillation  rate, 
there  are  also  times  when  scintillation  rate  is  reduced 
momentarily  and  the  average  amplitude  level  is  simulta¬ 
neously  decreased,  coinciding  with  the  increase  in  Faraday 
rotation  angle.  Two  such  cases  can  be  identified  in  Figure  5, 
at  0010  and  0038  LMT.  The  labeling  of  Faraday  rotation 
enhancement  cases  may  look  subjective,  especially  when 


they  occur  along  with  depletions  as  in  Figure  5.  Such  doubts 
may  be  further  reinforced  from  the  fact  that  the  in  situ  data 
have  clearly  shown  the  presence  of  isolated  depletions  in 
electron  density,  but  no  enhancement  cases  have  been 
reported  so  far.  However,  the  decrease  in  scintillation  rate 
and  scintillation  index  below  the  normal  level  at  the  time  of 
enhancement  and  the  rare  example  of  isolated  enhancement 
like  the  one  shown  in  Figure  6  make  us  believe  that  both 
depletions  and  enhancements  do  occur,  though  their  relative 
occurrence  may  be  affected  by  the  subjective  approach. 
Consequently,  the  typical  behavior  of  these  events  is  as 
follows;  The  depletions  in  TEC  are  usually  accompanied  by 
simultaneous  increases  in  fading  rate,  scintillation  index,  and 
the  average  amplitude  level,  while  the  enhancements  in  TEC 
are  accompanied  by  simultaneous  decreases  in  fading  rate, 
scintillation  index,  and  the  average  amplitude  level.  During 
the  1-year  period  from  September  1978  through  August  1979 
a  total  of  990  cases  consisting  of  710  isolated  depletions  and 
280  enhancements  have  been  identified.  Thus  the  cases  of 
enhancements  were  only  about  28%  of  the  total  cases.  The 
statistics  of  these  cases  follows. 

Figure  7  is  a  histogram  showing  the  occurrence  of  isolated 
depletions  and  enhancements  during  different  months  for  the 
period  September  1978  to  August  1979.  It  is  clear  that  the 
occurrence  of  both  depletions  and  enhancements  is  a  maxi¬ 
mum  in  the  month  of  October  and  a  minimum  in  August.  It  is 
noted  that  67%  of  depletions  occur  during  October-Decem- 
ber  and  80%  during  the  period  September-February,  while 
83%  of  enhancements  occur  during  September-December. 
Furthermore,  only  2%  of  depletions  and  4%  of  enhance¬ 
ments  occur  during  the  May-August  period.  Figure  8  shows 
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Fig.  7.  Histogram  showing  occurrence  of  isolated  depletions 
and  enhancements  during  different  months  for  the  l-year  period 
from  September  1978  to  August  1979. 


a  histogram  depicting  the  diurnal  variation  of  the  occurrence 
of  such  cases.  It  is  clear  that  these  depletions/enhancements 
occur  mostly  at  night,  between  1900  and  0300  LT,  with  a 
maximum  between  2000  and  2300  LT. 

A  histogram  depicting  the  duration  of  the  ionization  deple¬ 
tions/enhancements  is  given  in  Figure  9.  Most  of  these  cases 
last  less  than  20  min,  with  a  large  number  of  depletion  cases 
(I6S)  having  a  duration  between  8  and  10  min,  but  there  are 
also  a  number  of  depletions/enhancements  lasting  longer 
than  half  an  hour.  These  depletions/enhancements  have  a 
distribution  shown  in  Figure  10.  Most  of  the  cases  have  a 
decrease/increase  in  equivalent  vertical  electron  content  of 
less  than  4  x  I0'*  el/m^,  but  depletions  as  large  as  1.2  x  10'’ 
el/m’  have  also  been  observed,  which  is  about  1 1-17%  of  the 
total  diurnal  change  in  TEC  for  the  period  October-Decem- 
ber,  when  most  of  the  cases  occur. 

Modeling  of  the  Ionospheric  Bubble 

In  an  attempt  to  simulate  a  bubble  by  a  model  a  Gaussian 
lateral  distribution  of  electron  density  has  been  used.  As¬ 
suming  that  the  lateral  distribution  and  the  height  distribu- 
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Fig.  8.  Histogram  showing  the  diurnal  frequency  of  occurrence 
for  the  period  September  1978  to  August  1979. 


tion  enter  in  a  multiplicative  manner,  the  two-dimensional 
electron  density  can  be  expressed  as 

Nix,  h)  =  fih)  •  1 1  -  o  exp  -  ((.v  -  .v„)//|’|  ( 1 ) 

where 


h  height: 

X  lateral  distance; 
a  constant,  usually  smaller  than  unity; 

/  half  width  of  the  bubble; 

fih)  height  distribution  of  the  electron  density  of  the 
background  ionosphere. 

Equation  (I)  gives  an  electron  density  model  which  repre¬ 
sents  a  depletion  when  the  constant  a  is  positive  and  an 
enhancement  when  a  is  negative.  In  each  case  the  maximum 
deviation  occurs  when  x  =  xo.  When  ray  bending  is  ignored, 
the  optical  path  i/>  can  be  evaluated  by  the  integration 


j  n  dh  = 


>h» 


40.3 
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NihII  -  a  exp  -  |(.r  -  .r,i)//|’l 


C) 


where 

n  refractive  index; 
i/in  reference  phase; 

Njg  background  TEC.  equal  to  ffih)  dh. 

According  to  ray  theory  the  phase  variation  (2)  along  the 
wave  front  will  develop  into  an  amplitude  variation  after 
propagating  a  distance.  This  amplitude  variation  as  observed 
on  the  ground  can  be  shown  to  be  proportional  to  the  integral 
WehandUu.  1972] 
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Equation  (3)  gives  the  amplitude  variation  on  the  ground.  If 
the  amplitude  is  measured  in  decibels,  its  variation  is  then 
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Fig.  9.  Histogram  showing  the  duration  of  ionization  depletions 
and  enhancements  as  observed  along  a  fixed  path  for  the  period 
September  1978  to  August  1979. 


proportional  to  the  exponent  of  (3).  The  parameters  deduced 
from  Figure  4  and  used  to  calculate  an  amplitude  pattern  on 
the  ground  by  using  (3)  are  as  follows; 

rt  =  0.163  /=  136  MHz 

N,„  =  X  10'’  el/nr  /  =  IS.-S  km 

The  bottom  of  the  layer  was  assumed  to  be  at  a  height  of  300 
km.  Since  the  elevation  angle  is  28°,  the  path  S  turns  out  to 
be  640  km. 

The  calculated  amplitude  pattern  on  the  ground  is  depicted 
in  Figure  1 1 .  The  similarity  of  the  calculated  and  the 
experimentally  observed  average  amplitude  pattern  is  obvi¬ 
ous.  The  maximum  amplitude  is  dependent  on  the  width  of 
the  bubble  when  the  other  parameters  are  fixed.  In  our  case 
the  best  fit  was  achieved  at  21  =  37  km. 

The  minimum  points  of  amplitude  can  be  determined  quite 
accurately  from  the  observed  amplitude  variation  (Figure  4). 
The  time  difference  between  the  minima  of  the  observed 
amplitude  variation  is  approximately  450  s.  The  correspond¬ 
ing  length  in  the  calculated  amplitude  variation  (Figure  10)  is 
about  44  km,  giving  a  drift  velocity  of  about  98  m/s,  which 
agrees  well  with  drift  velocities  measured  independently  by 
other  means  (  Yeh  et  aL,  19811. 

As  was  pointed  our  previously,  the  amplitude  pattern  of 
enhancement  can  be  computed  simply  by  changing  the  sign 
of  a  in  ( 1).  On  a  logarithmic  scale  the  curve  will  be  the  .same 
shape  as  in  the  case  of  depletion,  but  turned  upside  down, 
because  in  (3).  a  is  in  the  exponent.  On  a  linear  scale  the 
nature  of  the  curve  will  be  different. 

It  is  obvious  that  the  suggested  Gaussian  model  is  not 
perfect.  However,  it  describes  quite  well  the  main  features 
of  the  average  amplitude  variation  caused  by  the  ionospheric 
bubble.  The  present  model  does  not  explain  the  change  in 
fading  rate  and  scintillation  index  caused  by  a  bubble  as 
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Fig.  10.  Histogram  showing  the  change  in  electron  content  for 
depletion  and  enhancement  cases  as  observed  during  the  peritHl 
September  1978  lo  August  1979. 
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Fig.  1 1 .  Amplitude  variation  of  a  radio  wave  on  the  ground  after 
passing  through  a  bubble.  The  width  of  the  bubble  is  37  km.  the 
frequency  136  MHz.  and  the  distance  from  the  layer  640  km.  The 
estimated  distance  between  the  two  minima  is  about  44  km. 


described  in  the  present  work.  This  is  because,  in  addition  to 
a  general  depletion  of  electron  density  inside  the  bubble  as 
modeled  in  ( I ),  there  are  also  sharp  fingerlike  wedges  which 
are  not  contained  in  (1).  Computations  by  Wernik  et  al. 
(1980]  show  that  regions  of  sharp  wedges  can  be  sources  of 
intense  scintillations.  The  model  constructed  by  them  is 
based  on  in  situ  measurements  [McClure  el  uL,  1977)  and 
shows  these  highly  complex  structures  with  steep  gradients 
in  the  center  of  the  bubble.  The  horizontal  gradients  gradual¬ 
ly  decrease  toward  the  edge  of  the  bubble,  as  also  revealed 
by  in  situ  measurements.  The  computed  pattern  indicates  an 
increase  of  the  scintillation  rate  with  increase  of  bubble 
thickness  and  shows  a  change  in  the  scintillation  rate  with 
increase  of  bubble  thickness  and  shows  a  change  in  the 
scintillation  rate  from  slow,  close  to  the  edges  of  the  bubble, 
to  fast,  at  its  center  |  VTerwA  eiai,  1980).  This  is  in  complete 
agreement  with  our  experimental  results. 
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